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Abstract 

We evaluate threshold resummed spectra in B ^ Xulv decays in next-to-leading order. We present results for 
the distribution in the hadronic variables Ex and m\ /E\ , for the distribution in Ex and for the distribution 
in Ex and Ei, where Ex and mx are the total energy and the invariant mass of the final hadronic state Xu 
respectively and Ei is the energy of the charged lepton. We explicitly show that all these spectra (where there 
is no integration over the hadronic energy) can be directly related to the photon spectrum in i? ^ Xg^ via 
short-distance coefficient functions. 
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1 Introduction and summary of the results 



A long-standing problem in particle physics is the understanding of strong interactions at low energies. While 
at very low energies, of the order of the hadronic scale A « 300 MeV, perturbative QCD is of no use and 
alternative methods have been developed in decades (such as quark models, chiral lagrangians, lattice QCD, 
etc.), at intermediate energies, of the order of a few GeV, perturbative computations can be combined with 
non-perturbative models to predict a variety of cross sections and decay rates. Among these moderate hard 
scale phenomena is beauty physics, which is indeed characterized by a hard scale of a few GeV. The measured 
decay spectra often receive large contributions at the endpoints — in the case of the hadron energy spectrum, 
in the middle of the domain — from long-distance effects related to soft interactions between the heavy quark 
and the light degrees of freedom. 

The main non perturbative effect is the well-known Fermi motion, which classically can be described as a 
small vibration of the heavy quark inside the B meson because of the momentum exchange with the valence 
quark; in the quantum theory it is also the virtuality of the heavy quark that matters. This effect is important 
in the end-point region, because it produces some smearing of the partonic spectra. 

These long distance effects manifest themselves in perturbation theory in the form of series of large infrared 
logarithms, coming from an "incomplete" cancellation of infrared divergencies in real and virtual diagrams. The 
probability for instance for a light quark produced in a process with a hard scale Q to evolve into a jet with an 
invariant mass smaller than m is written in leading order as 



= 1 - ^ "s log — 
2 \ 



Ai as 



1 da; dO^ 



(1) 



where co is the energy of a gluon emitted by the light quark normalized to the hard scale, 6 is its emission angle 
and Ai is a positive constant (see sec. OJ. The first integral on the r.h.s. is the real contribution while the 
second integral is the virtual one. Both integrals are separately divergent for w = — soft singularity — as 
well as for 9 = — coUinear singularity, but their sum is finite. "Complete" real-virtual cancellation occurs 
only for m = Q, i.e. in the completely inclusive evolution of the quark line, while for m < Q there is a left-over 
double logarithm because of the smaller integration region of the real diagrams. Multiple gluon emission occurs 
in higher orders of perturbation theory; it can be described as a classical branching process and gives rise to 
the double logarithmic series, i.e. to powers of the last term aglog^ [Q^/rn?^ on the r.h.s. of eq. JIE]. 

We may say that perturbation theory "signals" long-distance effects in a specific way — even though a 
quantitative description of the latter has to include also some truly non-perturbative component. A theoretical 
study of the universality of these long-distance effects can therefore be done inside perturbation theory, by 
comparing the logarithmic structure of different distributions. In other words, if these long-distance effects are 
universal, this has certainly to show up in perturbation theory: things have to work in perturbation theory first. 
The aim of this work is to study the relation of long-distance effects between different distributions by means 
of resummed perturbation theory. 

In general, let us consider the semi-inclusive decays 

B ^ Xq + (non QCD partons), (2) 

where Xq is any hadronic final state coming from the fragmentation of the light quark q = u,d,s and the 
non QCD partons are typically a photon, a lepton-neutrino pair, a lepton-antilepton pair, etc. This system 
of particle(s), with total four-momentum g^, constitutes a "probe" for the hadronic process, as in the case of 
deep-inelastic-scattering (DIS) of leptons off hadrons. Without any generality loss, we can work in the h rest 
frame, where p'^ = m.bv'^, with rub being the beauty mass and = (1; 0, 0, 0) being the classical 4- velocity. The 
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hadronic subprocess in ^ is characterized by the following three scales: 

rub, Ex and mx {mb > Ex), (3) 

where mx and Ex are the invariant mass and the total energy of the final hadronic state Xg, respectively. We 
are interested in the so-called threshold region, which can be defined in all generality as the one having 

mx <t: Ex. (4) 

The region (0J is sometimes called radiation-inhibited, because the emitted radiation naturally produces final 
states with an invariant mass of the order of the hard scale: mx ~ 0{Ex)- It is also called semi-inclusive 
because experimentally, to satisfy the constraint 10} , most hadronic final states have to be discarded. 

The processes we are going to consider are the well-known radiative decay with a real photon in the final 
state, 

B ^ X, +-f (5) 

and the semi-leptonic decay, 



B ^ Xu + I + ly. (7) 

In perturbative QCD, the hadronic subprocess in ^ consists of a heavy quark decaying into a light quark 
which evolves later into a jet of soft and coUinear partons because of infrared divergencies. In leading order, 
one only considers the emission of soft gluons at small angle by the light quark (see eq. (^); the final state Xq 
consists of a jet with the leading (i.e. most energetic) quark q originating the jet itself. In next-to-leading order 
one has to take into account two different single-logarithmic effects: (a) hard emission at small angle by the 
light quark q and (b) soft emission at large angle by the heavy quark. Because of (a), the final state consists 
of a jet with many hard partons and, in general, the leading parton is no longer the quark q which originated 
the jet itself. Because of (6), the final state does not contain only an isolated jet, but also soft partons in any 
space direction. The main result of is that the large threshold logarithms appearing in (O are conveniently 
organized as a series of the form: 

oo 2n ^2 
ra=lfc=l X 

= ci2 a{Q) log" + cn a(Q) log + C24 a'(Q) log^ + C23 a^iQ) log' ^ + , (8) 
mj^ mj^ mj^ mjr 

where a{Q) = as{Q) is the QCD coupling and the hard scale Q is determined by the final hadronic energy Ex 
5. 

Q = '2Ex. (9) 

These large logarithms are factorized into a universal QCD form factor. Let us summarize the derivation of 
lO and We take the infinite mass limit for the beauty quark while keeping the hadronic energy and the 
hadronic mass fixed ^: 

mb — !■ 00, with Ex and mx const. (10) 
"^The results for the semileptonic decay are easily extended to the radiative decay with the photon converting into a lepton pair, 

B ^ Xs + I + (6) 

^The factor two is inserted in such a way that the hard scale coincides with m;, in the radiative decay (see later). The essential 
point however is that Q is proportional to Ex via a proportionality constant of order one, whose precise value is irrelevant. 

®This limit has not to be confused with that one relevant for the shape function, also called structure function of the heavy 
flavors, which is Ex — * 00, mx — * 00 with m?^/Ex — > const (the latter implies m;, — > 00, but the converse is not true). The 
shape function describes soft interactions only and therefore does not factorizc the whole logarithmic structure, missing the large 
logarithms coming from hard coUinear emission off the light quark . 
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This takes us into an effective tlieory in which the beauty quarlc is replaced by a static quark, as recoil effects 
are neglected in the limit H1U|I . If we write the beauty quark momentum as pb = rribV + fc, where is a soft 
momentum, the infinite mass limit of the propagator is easily obtained as: 

where a = 'j^a'^. As discussed above, the beauty quark contributes to the QCD form factor via large logarithms 
coming from soft emissions, which are correctly described by a static quark. Since the light quark propagator 
is not touched by the limit we conclude that all soft and/or collinear emissions are correctly described 

by this limit. Since the heavy flavor mass has disappeared with the limit (|10|l . the only remaining scales in 
the hadronic subprocess are mx and Ex ■ Only one adimensional quantity can be constructed out of them, for 
example the ratio Ex/rrix, which is therefore the only possible argument for the large logarithms, in agreement 
with Furthermore, the hard scale Q is given by the greatest scale in the game, i.e. by the hadronic energy 
Ex, in agreement with 

The argument given above, however, is not rigorous: let us refine it. The limit (|1U|I is indeed singular in 
quantum field theory: one cannot remove degrees of freedom without paying some price. Let us consider for 
simplicity's sake the semileptonic decay O, even though the conclusions are general. The vector and axial- 
vector currents responsible for the 6 — > u transition are conserved or partially conserved in QCD, implying that 
the 0{a) virtual corrections are ultraviolet finite. These corrections contain however terms of the form 

70 a log ^ (ordinary QCD), (12) 

which diverge in the limit (|10|l (70 is a constant). If one takes the limit H1(J|) ab initio, i.e. before integrating 
the loop, some divergence is expected in the loop integrals, as it is indeed the case. Technically, that occurs 
because the static propagator is of the form 1/(^0 + ie) (see eq. (|11|) ') and, unlike the ordinary propagator, has no 
damping for |fc| — > cx). It can be shown that the h —> u vector and axial- vector currents are no more conserved 
or partially conserved in the static theory. Therefore, unlike the QCD case, the 0{a) virtual corrections are 
ultraviolet divergent in the static theory and produce, after renormalization, terms corresponding to (|12|l of the 
form 

7oQ;log-^ (effective theory), (13) 

Ex 

in which basically the heavy flavor mass m;, is replaced by the renormalization point — the coefficient 70 being 
the same. The hadronic subprocess in the static theory therefore has amplitudes depending on the physical 
scales Ex and mx as well as on the renormalization scale /i. If we neglect terms suppressed by inverse powers 
of the beauty mass ~ l/m^, we have that the physical scale rrifc is replaced by the renormalization point in 
the effective theory: rrib fi. The effective currents J,^ and the coupling constant a are renormalized at the 
scale /x: = Ji/(^) and a = a{^). The effective amplitudes contain terms of the form a" log'' m/^^ (k < n), 
which are large logarithms for fi » Ex or ^ <C Ex- To have convergence of the perturbative series, the large 
logarithms above must be resummed by taking fj. — 0{Ex), i-e. fj, = k Ex with k = 0(1). This implies that the 
effective currents and the coupling are evaluated at a scale of the order of the hadronic energy: = Ji,{k Ex) 
and a = a(k Ex)- We have therefore proved that the hard scale Q is fixed by the final hadronic energy Ex and 
not by the beauty mass m;,: 

Q = H = kEx with k = 0(1). (14) 



Let us go back to the general process Kinematics gives: 



2ii;x = m, ( 1 - 4 + ^ ) ■ (15) 

2 



The simplest processes are those with a light-like probe, i.e. with q = 0, where 



2Ex ^ mb[l + ^] ~ rub- (16) 

771 1, 
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This case corresponds to the radiative decay Q. In this case, the final hadronic energy is always large and of 
the order of the heavy-flavor mass: 



Q ^ nib 



(radiative decay). 



(17) 




(18) 



implying a substantial reduction of the hard scale: 



(19) 



This fact is one of the complications in the threshold resummation of the scmilcptonic decay spectra: while in 
the radiative decay the hard scale Q is always large in the threshold region, and of the order of mi,, this is 
no longer true in the semileptonic decay. The hadronic subprocesses have in general different hard scales in the 
two decays. If one integrates over q^, for example because of undetected neutrino momentum, there is a mixing 
of hadronic contributions with different hard scales in the semileptonic case. However, it turns out by explicit 
computation that the contributions from a large q^, i.e. with a small hard scale in the hadronic subprocess, are 
rather suppressed (see sec.0J. 

At fixed Q, the large logarithms in ^ can be factorized into a QCD form factor, which is universal in 
the sense that it depends only on the hadronic subprocess. The differences between, let us say, the radiative 
decay © and the semileptonic decay (0) only enter in the specific form of a short-distance coefficient function 
multiplying the QCD form factor (and in the form of a remainder function collecting non factorized, small 
contributions, see next section). 

The discussion above can be summarized as follows. The hard scale Q = 2 Ex in ^ appears in the argument 
in the infrared logarithms as well as in the argument of the running coupling. In the radiative decay, because of 
kinematics, the hard scale is always large and of the order of the beauty mass: Q ~ m?,, while in the semileptonic 
case kinematical configurations are possible with Q « m?, as well as with Q <^ m?,. The main complication 
in semileptonic decays is that by performing kinematical integrations (for example over the neutrino energy), 
one may integrate over the hard scale of the hadronic subprocess. While in radiative decays the hard scale is 
fixed, in the semileptonic decays there can be a mixing of different hadronic subprocesses. A non-trivial picture 
of some semileptonic decay spectra emerges: there are long-distance effects which cannot be extracted by the 
radiative decay, related to a small final hadronic energy, but their effect turns out to be small at the end because 
of a kinematical suppression of the states with a small hard scale. The decay spectra in (0 can therefore be 
divided into two classes: 

1. distributions in which the hadronic energy Ex is not integrated over. These distributions can be related via 
short-distance coefficients to the photon spectrum in the radiative decay (0). In particular, the structure 
of the threshold logarithms is the same as in decay ^ . In this paper we restrict ourselves to these simpler 
distributions; 

2. distributions in which the hadronic energy is integrated over and therefore all the hadronic energies 
contribute. These are for instance the hadron mass distribution or the charged lepton energy distribution. 
In all these cases, the structure of the threshold logarithms is different from that one in © and by far 
more complicated. The analysis of some of these distributions, which present novel features with respect 
to B ^ Xsj, is given in 0. 

Let us make a simple analogy with e^e^ annihilation into hadrons. In the center-of-mass (c.o.m.) frame, 
the final state consists of a qq pair, which are emitted back to back with a high virtuality and evolve later into 
two jets: 



e+ -I- e —> q + q Jq + Jq. 



(20) 



'^The same is also true for the radiative decay with the photon converting into a lepton pair JSJ. 
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Roughly speaking, the final state Xq in ||2J), consisting in a single jet, is "half" of that in H2UI) . consisting of 
the two jets Jq and Jq. Deviations from this independent fragmentation picture arise in next-to-leading order 
because of large- angle soft emission by the heavy quark in which has no analogue in l|2U|l . The structure 
of e^e~ hadronic final states is conveniently analyzed by means of so-called shape variables, one of the most 
studied being the heavy jet mass defined as 



max 



{ml,ml], (21) 



where m^j and to^ are the invariant masses of the particles in the right and left hemispheres of the event 
respectively. The hemispheres are defined cutting the space with a plane orthogonal to the thrust axis n, the 
latter defined as the direction maximizing 

El^^'-'^l' (22) 

i 

i.e. basically the sum of length of longitudinal momenta. The sum extends over all hadrons — partons in 
the perturbative computation. For mu <SC Q, where Q is the hard scale to be identified here with the c.o.m. 
energy, hard emission at large angle by the qq pair cannot occur and the final state consists of two narrow jets 
around the original qq direction, which can be identified with n. The 0{a) computation gives large logarithms 
of similar form to those in 0]: 

a(g)log2^ and a{Q) \og-\. (23) 
There is not a simple relation between, let us say, the heavy jet mass distribution at the Z'^ peak, 

'^'^ {Q = mz) (24) 



and the integral of this quantity over Q from a small energy e ^ niH up to mz with some weight function (j){Q): 

r dQm)^{Q). (25) 

Radiative B decays (jS)) and semileptonic spectra ((TJ in class 1. are the analog of the former distribution (|24|l . 
while semileptonic spectra in class 2. are the analog of the latter case (I25II . The analog of the suppression in 
the semileptonic spectra 2. of the contributions from large q^ is the suppression of the weight function 4'{Q) for 
Q < mz- 

Many properties of the distributions we are going to derive in this work can be understood with a qualitative 
discussion on the hadron energy spectrum, 

which exhibits a remarkable phenomenon related to the occurrence of infrared singularities inside the physical 
domain, instead than at the boundary as it is usually the case. This phenomenon has been studied in the 
framework of jet physics and is known as the "Sudakov shoulder" [HIEIEIEI- Let us discuss it in the present 
case in physical terms. In lowest order, the semileptonic decay {Tj) involves three massless partons in the final 
state: 

h ^ u + I + V. (27) 
According to kinematics, any final state parton can take at most half of the initial energy, implying that 

4°) =E^<^. (28) 

In lowest order, the final hadronic state consists indeed of the up quark only: X„ — u. To order a, a real gluon 
is radiated and the final hadronic state is a two-particle system: X„ = u -'r g. The final hadronic energy is not 
restricted anymore to half the beauty mass but can go up to the whole beauty mass: 

e'^x^ ^ Eu + Eg < rub. (29) 
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For example, just consider an energetic up quark recoiling against the gluon, with a soft electron and a soft 
neutrino. The relevant case for us is a final state with the up quark of energy « mb/2 and a soft and/or a 
coUinear gluon. Such a state has a total energy slightly above mb/2 and the matrix element is logarithmically 
enhanced because of the well-known infrared singularities. Such logarithmic enhancement cannot be cancelled 
by the 0(a) virtual corrections, because of their tree-level kinematical limitation H28|l. We are left therefore 
with large infrared logarithms, of the form 

alo,^{E.-^), alog(i^x-^) ^ ' (^0) 

which are final and produce an infrared divergence for Ex — > + mb/2. On the other hand, for Ex < 'ti(,/2 there 
are no large logarithms of the form alog'^(mf,/2 — Ex) {k — 1, 2), because in this case real-virtual cancellation 
may occur, and it actually does. Let us summarize: the 0{a) spectrum has an infrared singularity right in the 
middle of the domain, for Ex = Tnb/2, because the lowest order spectrum has a discontinuity in this point, 
above which it vanishes identically because of kinematics. 

This infrared singularity is integrable, as 

/ dEx a log"" [Ex < oo, (31) 

where 5 > is some energy-resolution parameter. The infrared divergence is therefore eliminated with some 
smearing over the hadronic energy, which experimentally is always the case. Furthermore, hadronization cor- 
rections certainly produce some smearing on the partonic final states because of parton recombination. In other 
words, non-perturbative mechanisms wash out this infrared divergence, which therefore does not present any 
problem of principle. As we are going to show, however, perturbation theory "saves itself" and no mechanism 
outside perturbation theory is needed to have a consistent prediction: resummation of the infrared logarithms 
in (|30(l to all orders completely eliminates the singularity, as in the cases of the usual infrared divergencies [H]. 
Since large logarithms occur for Ex ~ rub/i, we have that the hard scale is given for this spectrum by the 
beauty mass, 

Q = rrib (hadron energy spectrum), (32) 

just like in radiative decays. This equality is noticeable, as it comes from completely independent kinematics 
with respect to the one in (jSJ. There is therefore a pure short-distance relation between the hadron mass 
distribution in Q and the hadron energy distribution in Q. This property remains true when we consider 
non-perturbative Fermi motion effects, which are factorized by the well-known structure function of the heavy 
flavors, also called the shape function. 

This paper is organized as follows: 

In sec. ^ we presents the results for the resummed triple-differential distribution, which is the most general 
distribution and the starting point of our analysis; 

In sec. we review the theory of threshold resummation in heavy flavor decays, giving explicit formulas 
for the QCD form factor in next-to- next-to-leading order (NNLO). The transformation to Mellin space in order 
to solve the kinematical constraints for multiple soft emission is discussed, together with the inverse transform 
to the original momentum space; 

In sec. we derive the double distribution in the hadronic energy and in the ratio (hadronic mass) / (hadronic 
energy), which are the most convenient variables for threshold resummation (these are the variables w and u 
defined there) . The distribution in any hadronic variable can be obtained from this distribution by integration; 

In sec. (0 we present the results for the resummed hadron energy spectrum in next-to-leading order, whose 
main physical properties have already been anticipated here. We also compute the average hadronic energy to 
first order and compare with the radiative decay. The hadron energy spectrum with an upper cutoff on the 
hadron mass, which is the easiest thing to measure in experiments, is derived in leading order; 
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In sec. ^ we derive the double distribution in the hadron and electron energies, i.e. in the two independent 
energies. A peculiarity of this spectrum is that it is characterized by the presence of two different series of large 
logarithms, which are factorized by two different QCD form factors. Another peculiarity is that this double 
differential distribution contains partially-integrated QCD form factors instead of differential ones. That implies 
that the infrared singularities occurring in this distribution are integrable, as in the case of the Sudakov shoulder 
which we have discussed before; 

Finally, in sec. ^ we present our conclusions together with a discussion about natural developments. 



2 Triple differential distribution 



The triple differential distribution in the decay ((T)) is the starting point of our analysis. It has a resummed 
expression of the form 01 :® 

1 d^r 

— = C [x , w; a{w mi,)] a [u; a{wmb)] + d[x,u,w;a{wmb)] , (33) 

A. vAitJU\A> (AiiA/ yXJ 

where we have defined the fohowing kinematical variables: 

w = — - (0 < < 2), X = — ^ (0 < a; < 1) (34) 

mj, nib 



and ' 

(0 < w < 1). (35) 



Ex - ^/E\ - mj, 
Ex + V^x - 
It is convenient to express u as: 



with 



The inverse formula of (|36|l reads: 



(1+m)2 

The functions entering the r.h.s. of eq. H33() are: 



1 - ^/T^ 
1 + VI - 4y 



y = T^T (0 < 2/ < 1/4). (37) 



V = TT^- (38) 



• C\x,w\a(wm,h)\^ a short-distance, process-dependent coefficient function, whose explicit expression will 
be given later. It depends on two independent energies x and w and on the QCD coupling a; 

• a [m; Q;(it;TOfc)], a process- independent, long-distance dominated, QCD form factor. It factorizes the thresh- 
old logarithms appearing in the perturbative expansion. At order a: 

-(-;«) = ^(.)-^fi^) -^f-) +o(«^)' (39) 

*We have normalized the distribution to the radiativcly-corrected total semileptonic width V = 
To [l + qCf/tt (25/8 - 7r2/2) + 0{a'')] and not to the Born width PC'), as originally done in E|. We consider it to be a 
better choice because F, unlike r^'^\ is a physical quantity, directly measurable in the experiments and we are not interested in the 
prediction of total rates, but only in how a given rate distributes among different hadronic channels. 

^Note that a similar variable simplifies two-loop computations with heavy quarks 
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where Cp is the Casimir of the fundamental representation of SU{3)c, Cp — {N^ ^1)/ (2A^c) with Nc = 3 
(the number of colors) and the plus distributions are defined as usual as: 



P{u)+ = P{u)-5{u) / du'P{u'). 

'o 



The action on a test function f{u) is therefore: 



duP{u)+f{u)^ f duP{u)[f{u) ~ f{0)]. 

JQ 



(40) 



(41) 



The plus-distributions are sometimes called star-distributions and can also be defined as limits of ordinary 
functions as: 



P{uh 



lim 



= lim 

£-.0+ 



(u- e)P(u) - 5{u) / du' P{u') 



e)P{u) ^ 5(u- e) / du P{u) 



lim — 

£-►0+ du 



u-e) I du' P(u') 



(42) 



An important property of the plus-distributions is that their integral on the unit interval vanishes: 



[ P{u)+ du = 0. 
Jo 



(43) 



We have assumed a minimal factorization scheme in eq. (|39|l . in which only terms containing plus- 
distributions are included in the form factor. The resummation of the logarithmically enhanced terms in 
a to all orders in perturbation theory will be discussed in the next section; 

d{x,u,w; a) is a short-distance, process-dependent, remainder function, not containing large logarithms. 
Formally, it can have at most an integrable singularity for u — > -1-0, i.e. we require that: 



lim / du'd(x,w,u':a) = 0. 



(44) 



This term is added to C • cr in order to correctly describe the region u ~ 0(1) and to reproduce the total 
rate. It depends on all the kinematical variables x, w and u and the explicit expression will be given later. 



Eq. (|33|l is a generalization of the threshold resummation formula for the radiative decay in Q [51 11U|: 

1 dTn 



Tr dts 



Cfl[Q;(u;TOb)] cT[ts] a{wmb)\ + [t^; a{wmb)] ~ 



where^° 



m 



(45) 



(46) 



In this simpler case, the coefficient function Cji{a) does not depend on any kinematical variable but only on 
the QCD coupling a and has an expansion of the form^^: 



(47) 



^"The relation with the photon energy x-f = 2_E-.y/r?i(, is ts = 1 — x^. 

^^We perform expansions in powers of a, while the traditional expansion is in powers of a/(27r). 
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where C^^ are numerical coefBcients. Basically, going from the 2-body decay Q to the 3-body decay Q, the 
coefficient function acquires a dependence on the additional kincmatical variables, namely two energies. The 
remainder function in eq. 145(1 depends on the (unique) variable tg and has an expansion of the form: 



2^(2), 



(48) 



The main point is that the QCD form factor a in the same in both distributions and H45|l. explicitly showing 
universality of long-distance effects in the two different decays. By universality we mean that we have the same 
function, evaluated at the argument u in the semileptonic case and at tg in the radiative decay. This property 
is not explicit in the original formulation "11", in which the form factors differ in subleading order (see next 
section). 

Since, as shown in the introduction, w 1 in the radiative decay, we can make everywhere in eq. (|45|) the 
replacement 

a{w m}j) — !■ a{mi,) (radiative case only), (49) 



to obtain: 



1 dTn, 
Vr dts 



Cr [a{mb)] a [tg] a{mb)] + dR [tg; a{mb)] . 



(50) 



The distribution contains now a constant coupling, independent on the kinematics a{'mh) ~ 0.22. The replace- 
ment (|49|) cannot be done in the semileptonic case. 

In |2j the triple differential distribution was originally given in terms of the variable y instead of w, with the 
latter it = 1 — ^ being introduced in The variables u and y coincide in the threshold region in leading twist, 
i.e. at leading order in u in the expansion for u — > 0, as y — u + 0{u^). Going from the variable y to the 
variable u only modifies the remainder function. The advantages of u over y are both technical and physical: 

• u has, unlike y, unitary range; 

• when we impose the kinematical relation between hadronic energy Ex^ and hadronic mass mx, of the 
radiative decay lO, u exactly equals tg: 



l-\Ex,=mb/2{l+mj, /ml) — 



(51) 



This property suggests that some higher-twist effects may cancel in taking proper ratios of radiative and 
semileptonic spectra. 



Let us now give the explicit expression of the coefficient function in the semileptonic case: 
where 



(52) 



C(i)(x, w) 



I2{w - x){l +x-~w); 

12—{w-x)l (l + x-w) 
I 

X log w \ 
+ 2(T^f 



T • / N , 1 /-, \ 3 , wlogw 35 

Li2 u) + log u; log 1 ~w)~ - \ogw - — - — 

2 2(1 — w) 8 



(53) 



(54) 



with X — \ — X ^"^ . Note that the coefficient function contains the overall factor w ~ x — x^, which vanishes 
linearly at the endpoint of the neutrino spectrum. We have defined Xi, ^ I — Xi, and Xi, = 2Ep /nib. 

^^To avoid spurious imaginary parts for w > 1 one can use the relation Li2(ji') = — Li2(l — w) — logu)log(l — ui) + 7r^/6. 
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Unlike the coefficient function, the remainder function d{x,u,w,a) has an expansion starting at 0{a): 

d (x, w, u; a) = a d'^^'' (x, w, u) + rf*-^-* {x, w, u) + 0{a'^). 
Omitting the overall factor Cf/tt, we obtain: 



(55) 



d^^\xTWTu) = 



+ 



+ 



-Sw'^ {2A + 3w-8x) 9w^ {24: + 3w~8x) 



4(1 +uy 



8(l + < 



9 (-12 + w) (-2 + wf {w -2xf 9 (-12 + w) (-2 + wf [w - 2xf 



16 (1 - uf 



32(1 -u)' 



iw'^ (32 - 47w - 8m;2 + 16 x + 20 w x + ^ + 8 -iwx'^) 

iw'^ (-64 + 94w + 40w2 + 3u,3 _ 32x - 40w;x - lOw^x - IQx"^ + Qwx'^) 

8(l + w)^ 



+ 



(640 w - 368 - 200 - 16 w'' + 3 - 384 x + i2Qwx + 528 x 

1 + u V 



64 (1 + u) 

112 w^x- 16w''x- 256^2 - 48wx^ - 224 w^^^ ^ 2Aw^x'^ 
3 



256 w + 528 - 200 w'' - 16w^ + Sw" + 512 x - 1472 wx 



64 (1 - u) 

+ 528w^x + 112w^x- 16w^x + 640x^ - 48wx^ - 224:w'^x'^ + 24w^x'^ 

9fi;^logu Qui^logu 9 (— 12 + (— 2 + w)^ (ui — 2 x)^ log w 
A{l + uf ^ 2(l + u)^ 16(1 -u)^ ^ 

9 (-12 + w) (-2 + wf {w - 2xf log u 
16 (1 - uf 

(-10+16w + w2 + 8x-2wx-2x2) log u 



8(1 + uY 

(-10+ 16u; + 7w;2 + 8x- 2wx- 2^2) log 



3 log u 



TT w ( - 144w + 208w^ + 16w^ + w"* - 64x- 80wx - IGw^x^ 

64(l + u)^ V 

Sw^x + ASx'^ -mwx'^ + 16w^x^'j + 

^ " , f- 256 w + 624 -304^3 + 16^4 + + 512^- 1856 tux 
64(1 -uy V 

+ 944 X - 16 X - 8 w"* X + 1024 - 464 wx^ - 96 x^ + 16 x^^ 



(57) 



The remainder function is a combination of rational functions of u multiplied in some cases by logu, with 
coefficients given by polynomials in w and x. 

^^The 0{a) function is obtained from that one given in |2] dPi\ {^,'W,y) in terms of the variables z = 1 — y and f = 1 — 4y, by 
using a relation extending eq. (23) of IIUI : 



d'l' {x,w,u) = d[]l{x,w,y{u)) ^{u) + C(°'(x,ii)) 



Cf 



log M + 7/4 \ogy{u) + 7/4 dy 
u yiu) du 



(u) 



(56) 



This function can also be obtained with a direct matching with the 0{a) triple differential distribution computed in ^ after a 
change of variable (see the end of this section for a discussion about matching) . 
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The main point about the semileptonic decay {T)) is that it has — unhke the radiative decay ||SJ) — 7^ 
and consequently the form factor depends not only on u but also on the hadronic energy w through the running 
coupling: 

a = a[u; a{wmb)]. (58) 
The form factor is therefore a function of two variables. 

We work in next-to-leading order (NLO), in which only the 0(a) corrections to the coefficient function and 
remainder function are retained (see next section). Since the difference between a(wnii,) and a{mi,) is 0{a^), 
we can set w = 1 in the argument of the coupling entering the coefficient function and the remainder function. 
We then obtain the simpler expression: 

1 <fr 

C [x,w; 0(7711,)] a [u; a{wmb)] + d[x,u,w; a{mb)] (NLO). (59) 



r dxdudw 



Note that we cannot set w = 1 in the coupling entering the form factor, because in the latter case a is multiplied 
by large logarithms, which "amplify" 0{a^) differences in the couplings (see next section). 

Let us make a few remarks about the final result of this section, eq. (|59|l : 



• it describes semi-inclusive decays, in which the internal structure of the hadronic final states is not ob- 
served, but only the total mass and energy are measured. Less inclusive quantities, such as for instance 
the energy distribution of the final up quark (i.e. the fragmentation function of the up quark), cannot be 
computed in this framework; 

• it constitutes an improvement of the fixed-order 0(a) result in all the cases in which there are large 
threshold logarithms. In all the other cases, where there are no threshold logarithms, such as for example 
the dilepton mass distribution ^12| , there is not any advantage of the resummed formula over the fixed- 
order one. 



In the next sections we integrate the resummed triple-differential distribution to obtain double and single 
(resummed) spectra. There are two methods to accomplish this task which are completely equivalent: 



1. The first method involves the direct integration of the complete triple-differential distribution. Schemati- 
cally: 

(spectrum) = C -a + Id. (60) 



Large logarithms come only from the first term on the r.h.s. of H60|l . while non-logarithmic, "small" terms 
come both from the first and the second term. To obtain a factorized form for the spectrum analogous 
to the one for the triple-distribution, in which the remainder function collects all the small terms, one 
rearranges the r.h.s. of (|6()|l : the small terms coming from the integration of O • cr are put in the remainder 
function; 

In the second method, one integrates the block C ■ a only and drops the small terms coming from the 
integration. The remainder function is obtained by expanding the resummed expression in powers of a 
and comparing with the fixed-order spectrum. 



3 Threshold Resummation 

It is convenient to define the partially integrated or cumulative form factor I](m, a): 



S(u; a) — I du a{u ; a). (61) 
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Performing the integrations, one obtains for the 0{a) form factor: 



S(K;a) = l-^i^ + ^L + 0(a^), (62) 

2 TT 4 TT 

where 

L = log-. (63) 
u 

E contains a double logarithm coming from the overlap of the soft and the coUinear region and a single logarithm 
of soft or coUinear origin. The normalization condition reads: 

E(l; q) = f dua{u;a) ^ I. (64) 

"'0 

As already noted, we have assumed a minimal factorization scheme, in which only logarithms and not constants 
or other functions are contained in the form factor. The expression of the partially integrated form factor E is 
technically simpler than the one for the differential form factor a, as it involves ordinary functions instead of 
generalized ones. Furthermore, in experiments one always measures some integral of a around a central u value 
because of the binning. 

In the limit w — > 0+, no final states are included in the distribution and therefore one expects, on physical 
grounds, that 

lim T,{u; a) ^ 0. (65) 

The 0(a) expression (|62|l does not have this limit and it is actually divergent to — oo — a completely un-physical 
result. In general, a truncated expansion in powers of a is divergent for m — > 0"*", because the coefficients diverge 
in this limit. Therefore, one has to resum the infrared logarithms, i.e. the terms of the form a" L'^, to all orders 
in perturbation theory. In higher orders, S contains at most two logarithms for each power of a, one of soft 
origin and another one of coUinear origin. Its general expression is then: 

Y.{L,a) = \ + (66) 

Tl=l k = l 

where Y^nk are numerical coefficients. At present, a complete resummation of all the logarithmically-enhanced 
terms on the r.h.s. of eq. (|66|) is not feasible in QCD: one has to resort to approximate schemes. The most 
crude approximation consists of picking up the most singular term for w — > 0+ for each power of a, i.e. all the 
terms of the form: 

a" L^" (double logarithmic approximation). (67) 

In this approximation, we can neglect running coupling effects and effects related to the kinematical constraints: 
higher orders simply exponentiate the 0(a) double logarithm and one obtains 

E (u; a) — e~ "Z'^'^) ^ (double logarithmic approximation). (68) 

Let us note that the resummed expression unhke the fixed-order one does satisfy the condition (|^ . 
The exponent in the resummed form factor involves a single term, — Cp aj (27r) L^, and has therefore a simpler 
form than the form factor itself. This remains true when more accurate resummation schemes are constructed, 
so it is convenient to define G as: 

E = e^. (69) 
It can be shown that the expansion for the function G is of the form |13|: 

oc n+1 
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where Gnk a-re numerical coefficients. Let us note tliat tlie sum over k extends up to n + 1 in (|70|) . while it 
extends up to 2n in the form factor in eq. This property is a generalization of the simple exponentiation 

of the 0{a) logarithms which holds in QED and is called generalized exponentiation. In general, this property 
holds for quantities analogous to the semi-inclusive form factors, in which the gluon radiation is not directly 
observed. One sums therefore over all possible final states coming from the evolution of the emitted gluons 
(inclusive gluon decay quantities). The property expressed by eq. (|70|l does not hold for quantities in which 
gluon radiation is observed directly, as for example in parton multiplicities, where different evolutions of gluon 
jets give rise to different multiplicities. 

3.1 A^-space 

A systematic resummation is consistently done in TV-moment space or Mellin space, in which kinematical 
constraints are factorized in the soft limit and are easily integrated over |14j . One considers the Mellin transform 
of the form factor (t(u; a): 

aN{a) = / rfu(l-it)^"V(u; a). (71) 

The threshold region is studied in moment space by taking the limit N ^ oo, because for large N the integral 
above takes contributions mainly from the region u ^ 1. For example, the Mellin transform of the spectrum in 
eq. H5()(l is of the form 

(1 - t, f-' — —^dU = Cnia) a^ia) + ^^.^(a), (72) 
I R ats 

where 

dR,Nia) ^0 for iV ^ c». (73) 

The total rate in Mellin space is obtained by taking iV = 1. 

It can be shown ^1 ^ E| that the form factor in iV-space has the following exponential structure: 

(Tat (a) = e<="("\ (74) 



where 



z 



1- z 



G^{a) = I dz—^ — { I ^ ^ [a(fc2)] -|- B [a{Q^{l - z))] + D [a{Q\l - zf)] ) . (75) 



'(1-2) 



Let us note that a prescription has to be assigned to this formula since it involves integrations over the Lan- 
dau pole The functions entering the resummation formula have a standard fixed-order expansion, with 
numerical coefficients: 

oo 

A{a) = ^ A„ a" = Aia + A^o? -I- A-^o? + A^a^ + ■■■ (76) 

n=l 

oc 

B{a) = '^B^a'' ^ Bia + B2a^ + Bsa^ + ■ ■ ■ (77) 

n=l 

oo 

D{a) = Dn a'^ = Dia + D2a^ + Dsa^ + ■ ■ ■ (78) 
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The known values for the resummation constants read: 



Bi 
B2 

D2 



1 

Cf_ 
2 



Ca 



p2 



67 z{2) 
36 

245 
"96" 



2 
11 



18 



(79) 
(80) 



^(3) - — z(2) + — z(A)] - CAUfi— + — z(3) -—z(2)] + 
24 ^ ^ 36 ^ ^ 8 ^ V ^432 12 ^ 18 ^ V 



/55 z(3)x ^2 
^^"H96'~J" 108 



_3Cf 

4 TT 

^2 



Ca 



3155 11 5 

h — z(2 + - z 3 

864 12 ^ 2 ^ ' 



Cf 



Cf 



io8-r(^)+ — 



1L 
54 



3 3 3 /247 z 2 

— + -z 3 - - z(2 + rif — 

32 2 ^ ^ 4 ^ V V 432 6 



(81) 

(82) 
;(83) 
(84) 
(85) 



where Ca — Nc = 3 is the Casimh of the adjoint representation. The coefficients Ai, Bi and Di are renormaliza- 
tion-scheme independent, as they can be obtained from tree-level amplitudes with one-gluon emission (see later). 
The higher-order coefficients are instead renormalization-scheme dependent and are given in the MS scheme 
for the coupling constant 

To this approximation, the first three orders of the /3-function are also needed [T^l^ : 





1 






An 






1 








247r2 


17 


f32 = 


1 




r2857 


647r3 


[ 54 



5CA + 3CF)nf 



(86) 
(87) 



/1415 205 



79 



11 



C^]nf+[f^CA + ^CF]nj 
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As is well known, f3a and /3i are renormalization-scheme independent, while P2 is not and has been given in the 
MS scheme. We define the /3-function with an overall minus sign: 



da 



-/3(a) 



/3o a - /3i - (32 ct 



d log ^'^ 

The running coupling reads: 

1 /3i log(log/iVA') , /32 1og' (log/xVA^) - log(logAiVA2) - 1 



/3ologAi2/A2 /33 logVVA^ 



log^ ^2/^2 



(89) 



f32 



/3oMogVVA2 



The functions A{a), B{a) and D{a) have the following physical interpretation (see for example |21ll22j '): 



(90) 



The function A{a) involves a double integration over the transverse momentum kt and the energy w of 
the emitted gluon and represents emissions at small angle and at small energy from the light quark. The 



^''A discussion about the scheme dependence of the higher order coefficients A2, B2, etc. on the coupling constant can be found 

1 El. 
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leading term Ai is the coefficient of tliat piece of the matrix element squared for one real gluon emission, 
which is singular in the small angle and small energy limit: 



duj dO'^ , duj dk? 
A,a-^^A,a-^, (91) 



UJ 



where kt ~ w is the transverse momentum of the gluon. In H91|l we have given the representation of the 
integral both in the angle 9 and in the transverse momentum kt- The subleading coefficients A2, A3, etc. 
represent corrections to the basic double-logarithmic emission. The function A{a) "counts" the number 
of light quark jets in different processes, i.e. we can write 

= nqA{a), (92) 

where Uq is the number of primary light quarks in the process P. For example, in e+e^ annihilation 
into hadrons Uq — 2, while in the heavy flavor decays ^ Uq — \. Since soft gluons only couple to the 
four- momentum of their emitters and not to their spin, the function Ag{a) for gluon jets is obtained from 
the quark one A{a) simply taking into account the change in the color charge, i.e. multiplying by Ca/Cp 
123; 

• the function B{a) represents emissions at small angle with a large energy from the light quark. Bi is the 
coefficient of that piece of the matrix element squared which is singular in the small angle limit: 

, , dk'l , , 

Biadio—r^ ^ Biaduj —Ty-. (93) 
6'' kf 

The non logarithmic integration over the gluon energy lu has been done and does not appear explicitly in 
eq. (|75|) : the integration over the angle or the transverse momentum kt is rewritten as an integral over 
z. The function B{a) counts the number of final-quark jets, i.e. 

B'^P\a) = niB{a), (94) 

where n; is the number of primary final quarks in the process P. For example in e~^e~ annihilation into 
hadrons = 2, while in DIS or in the heavy flavor decays ^ ni = I. Since hard coUinear emissions are 
sensitive to the spin of the emitting particles, the gluon function Bg{a) is not simply related to the quark 
one B{a) 

• the function D{a) represents emissions at large angle and small energy from the heavy quark. Di is the 
coefficient of that piece of the matrix element squared which is singular in the small energy limit: 

Dia — d0\ (95) 
to 

The non logarithmic integration over the angle or the transverse momentum kt has been done and does 
not appear explicitly in eq. (|75|l : the integration over the energy lu is rewritten as an integral over z. 
Di = in all the processes involving light partons only, as for instance DIS, Drell-Yan (DY) or e+e" 
annihilation into hadrons, while it is not zero in all the processes containing at least one heavy quark, 
such as for example the heavy flavor decays (0) . Note that the effective coupling appearing in the D terms 
is a [Q^{1 — z)^] and is therefore substantially larger for 1 — z ^ 1 than the coupling entering the hard 
collinear terms, namely a [Q^(l — z)] . 



Eq. (|75|l is therefore a generalization of the 0{a) result, possessing a double logarithm coming from the overlap 
of the soft and the collinear region and a single logarithm of soft or collinear origin (see eqs. 1)91(1 . I(93|l and 
195(1 ') The functions A{a) and -6(0;) are believed to by universal, i.e. process independent to any order in 
perturbation theory, as they represent the development of a parton into a jet, i.e. one-particle properties. The 
function £'(0;) on the contrary is process-dependent, as it describes soft emission at large angle, with interference 

^^Let us remember however that only two of the three functions appearing in eq. 1751 are independent llbl . 



15 



contributions from all the hard partons in the process, i.e. it describes global properties of the hadronic final 
states. Let us observe that A2 and D2, unlike B2, do not have a Cp contribution. That is a consequence of the 
eikonal identity, which holds in the soft limit According to this identity, the abelian contributions simply 
exponentiate the lowest order 0(a Cp) term, just like in QED. That means that there are no higher order terms 
in the exponent Gn- Because of similar reasonings, A3 does not have a Cp contribution. 

Despite its supposed asymptotic nature, the numerical values of the coefficients show a rather good conver- 
gence of the perturbative series. Note that all the double-logarithmic coefficients Ai are positive, implying an 
increasing suppression with the order of the expansion (up to the third one) of the rate in the threshold region. 



On the contrary, the single-logarithmic coefficients Bi and Di — with the exception of B2 — are all negative 
and therefore tend to enhance the rate in the threshold region [23 . We have: 

Ai = +0.424413; (96) 

A2 = +0.420947 - 0.0375264 n/ = 0.308367; (97) 

A3 = +0.592067 - 0.0923137n/ - 0.000398167 = 0.311542; (98) 

Bi = -0.318310; (99) 

B2 ^ +0.229655 + 0.04020 = 0.350269; (100) 

Di = -0.424413; (101) 

D2 = -0.556416 + 0.002502 n/ = -0.548911. (102) 

With our definition, the /3-function coefficients are, as well known, all positive. 

/3o = +0.87535 - 0.05305 n/ = +0.71620; (103) 

Pi = +0.64592 - 0.08021 n/ = +0.40529; (104) 

P2 = +0.71986 - 0.140904 n/ + 0.003032 = +0.324436. (105) 



In the last member we have assumed 3 active flavors (ny = 3). 

Let us now discuss the computation of the coefficients entering the resummation formula. The occurrence of a 
Sudakov form factor in semileptonic B decays was acknowledged originally in j2Sl, where a simple exponentiation 
involving Ai and Bi -\- Di was performed. The coefficient A2 was computed for the first time, as far as we know, 
in j26| . It was denoted Ai K since it was considered a kind of rcnormalization of the lowest-order contribution: 

Aia ^ Aiail + K a). (106) 

The coefficient A2 was obtained from the soft-singular part of the q ~* q two- loop splitting function |2I| , that 
is as the coefficient of the 1/(1 — z) term^^. A2 was subsequently recomputed in in the framework of 
Wilson line theory, where the function A{a) has a geometrical meaning: it is the anomalous dimension of a 
cusp operator, representing the radiation emitted because of a sudden change of velocity of a heavy quark, 

oc 

rctisp(<^) = r^„gp a = r^Jgp a + r^J^po; + ••• (lo^) 

n=l 

Indeed, it has been explicitly checked up to second order that these two functions coincide: 

A{a) = T,usp{a)- (108) 

Let us note that: 

• the theory of Wilson lines and Wilson loops; 

^^This is exactly the same procedure which has been followed to derive the third-order coefficient A3 from the three-loop splitting 
function 
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• the eikonal or soft approximation in perturbative QCD; 

• the heavy quark effective theory (HQET) and the large energy effective theory (LEET), 



aU involve basically the same structure, i.e. the same propagators and vertices and the same amplitudes. Since 
the same structure has been studied in different frameworks, there is multiple notation and terminology for the 
same objects. Let us stress however that in ordinary QCD the function A{a) is not an anomalous dimension, 
since it is not obtained from ultraviolet 1 /e poles in renormalization constants but from infrared poles or from 
finite parts of scattering amplitudes. A{a) becomes an anomalous dimension in the effective theory because 
the latter has additional ultraviolet divergencies with respect to QCD. While in QCD one has to subtract only 
ultraviolet divergencies related to coupling constant renormalization, in the effective theory one has also to 
subtract additional ultraviolet divergencies related to the cusp operators. A scheme dependence is therefore 
introduced in the effective theory, which is not present in full QCD. It seems to us therefore that the equality 
I|1U8|) is not guaranteed a priori in higher orders and may require a specific scheme for the subtractions in the 
effective theory. At present, A3 has only been derived in full QCD and not in the effective theory. 

The coefficient B2 has been computed by means of the second order correction to the inclusive DIS cross 
section, which contains the combination B2 + DJ^^^ (the DIS analogue of eq. \i'A<6\ , see later) and by means of 
the third order correction, which contains the different combination B2 + 2D2^^ (the DIS analogue of eq. H139() , 
see later). The knowledge of the fermionic contribution to the 0{a^) DIS cross section was sufficient for a 
complete determination of B2 •''>() , with later checks offered by the complete computation |28l IM] . 

An incorrect value for the coefficient D2 for heavy favor decays has been obtained in the original computation 
in 1221 J where the technique to compute real and virtual diagrams in the effective-theory in configuration space 
has also been developed. In the coefficient of the single logarithm in the radiative decay © to order has 
been presented, from which the correct value of D2 can be extracted (let us note however that numerically the 
two values are not very different). A second order computation of heavy flavor fragmentation in ordinary QCD 
was presented in [SJ, which allows the determination of the sum B2 + D^^""^ (the analogue of eq. (|136|l . see 
later). Using an identity relating the coefficient for heavy flavor fragmentation with that one for heavy flavor 
decays, and subtracting the known value for the universal coefficient B2^ the correct value for D2 was explicitly 
derived in (see also jSHl)- Still in by repeating the Wilson line computation of errors were found 
and the same value of D2 extracted from heavy flavor fragmentation was re-obtained. Recently, the second 
order contribution of the chromomagnetic operator O7 to the photon spectrum in the radiative decay ((SJ has 
been calculated 37 , conffiming these results (see also j38j). 

According to the previous remarks concerning the relation between A{a) and Tcusp{oi), we believe it is a 
non-trivial fact that the same value of D2 is obtained with two completely different methods: 

• a direct computation in the effective theory, which describes the soft region only; 

• an extraction from an ordinary QCD computation, which gives the sum of the soft and the coUinear 
contributions B2 + D2, by subtracting the coUinear contribution B2 obtained from second order and third 
order DIS computations. 

The following expansion holds true for the exponent: 



The expansion of the logarithm of the form factor has a similar structure in physical space and in A^-space; 
roughly speaking, going to A^-space, logl/u logN. 




(109) 



where 



I = logN. 



(110) 
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As already discussed, we are interested in the large- TV limit; the 0{1/N) terms can be neglected in our 
leading-twist analysis. A resummation of all the logarithmically-enhanced terms in (|109|l is at present unfeasible 
in QCD even in iV-space, so one has to rely on approximate schemes. Let us discuss the fixed-logarithmic 
accuracy scheme: 

• Leading order (LO). One keeps in the exponent Gat (a) only the leading power of the logarithm for each 
power of a, i.e. fc = n + 1: 

oo 

G^o ^ ^G„„+ia"r+i = Gi2af + G23 a' + O(a'). (HI) 

n=l 

The coefficient function is kept in lowest order, i.e. G^*^ = 1 and the remainder function is completely 
neglected, i.e. d^'~' = 0; 

• Next-to-leading order (NLO). One keeps in Gat (a) also the terms with n — k, i.e.: 

00 

Gn"-" = ^ [G„„+ia"P+i + G„„a"P] 

n=l 

= Gi2af + Glial + G23a^l^ + G22a'^l^ + 0{a^). (112) 



To 0(a) one retains both the double and the single logarithm. In general for each order in a one keeps 
the principal two logarithms. One also keeps the 0{a) terms both in the coefficient function and in the 
remainder function: 

G^^o = l + aG«; d^^^ = ad«. (113) 

The one-loop coefficient function is needed because of the factorized form of the QCD form factor. One 
has indeed a resummed expression of the form: 



1 + aGW 



By expanding the exponent in powers of a, a term coupling the coefficient function and the double 
logarithm is obtained: 

a'^G^^^Guf, (115) 
which must be included in the NLO approximation; 

• Next-to-next-to-leading order (NNLO). One keeps in Gat also the terms with n = fc — 1, i.e.: 

00 

N = 2^[Gnn+ia I +Gnna I +Gnn-ia I \ 

n=\ 

= Gi2af + Glial + G23aH^ + Gaaa'/' + Gaia'/ + 034^^/* + G^saH^ + Ggaa^^' + 
+ 0{a'^). (116) 

To 0{a^), all the infrared logarithms are included. In general, for each order in a, one keeps the principal 
three logarithms. The first omitted term is the single logarithm to order a^. One has also to keep the 
0{a'^) terms both in the coefficient function and in the remainder function: 

G^^^o = 1 + aG(i) + d^^^o ^ ^ ^(i) ^ ^2 ^(2) ^ (^^7) 

The classes of logarithms discussed above can be explicitly resummed by means of a function series expansion 
of GAr(a) P: 

00 

GNia) = IgiiX) + 

a"g2+n(A) = I gi{X) + g2{\) + a gsiX) + a"^ g4{X) ~\ , (118) 
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where 

A ^ Poal. (119) 
The gi{X) are homogeneous functions of A and have a series expansion around A = 0: 



g^{X) = Y.9^n^"■ (120) 



n=l 

In LO one needs only the function gi, in NLO one need also g2, in NNLO also 53 is needed and so on. The 
explicit expressions read: 

51(A) ^ _-^[(i_2A)log(l-2A)-2 (l-A)log(l-A)]; (121) 



A, 13, 



52(A) = log(l - 2 A) + log(l - A) + — ^ [log(l - 2 A) - 2 log(l - A)] + (122) 

^ Po Po 2 po 



^ [2 1og(l-2A) + log2(l-2A)-4 1og(l-A)-21og2(l-A)] + 
4po 

+ ^ [log(l - 2 A) - log(l - A)] + [log(l - 2 A) - 2 log(l - A)] log ^. 
Po 2 /i^o 

The function 51(A) in jS] is in agreement with that one obtained originally in ^Jj. 52(A) in [3] differs instead from 
the corresponding 5|'(A) obtained in and it is equal to the corresponding function entering the B Xs^ 
spectrum; the formalism we use makes explicit the universality of soft gluon dynamics in semileptonic and 
radiative decays. The NNLO function 53 has the rather lengthy expression: 



ff3(A) 



Da A 2Di"fEX Dif3i f 2X log(l-2A)\ B2X Bije X 



/3o(l-2A) 1-2A 2/3o^ V1-2A ' 1-2A ; /3o (1 - A) 1-A 

B^g ( X log(l-A) \ A3 / A A \ ( 1 l_ 

/^o^^Hl-A 1-A ) 2/3o' V1-2A l-x) /3o \l-2X 1-A 

A2 /3i / 3 A 3 A log(l - 2 A) 2 log(l - A) 



2/3o^Vl-2A 1-A' 1-2A 1-A 
Ai7b2 / 4A A A v4i7r2 / 4A A 



2 V1-2A 1-Ay 12 V1-2A 1-A 
Ai/3i7B / 1 1 , log(l-2A) log(l-A) 



/3o2 \^1-2A 1-A' 1-2A 1-A 
DiX 11'^ BiX ji'^ A2 ( X X \ jp- 



/?o (1-2A) ^°^Q2 (1 „ A) ^^2(^i_2A l-Aj^°^Q2 

Ai7g / 2A A \ 

A U-2A l-Ay* °®Q2 + 
AiA ^ A A log(l-2A) , Alog(l-2A) Alog(l-A) ^ , 

^ (vi^^T^ 2 + 1-2A T^^j 

^^^^'^ '^Vog^^. (123) 



2/3o V1-2A I-A; Q2- 
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The function (73(A) was originally computed in |39| . where the first NNLO resummation in heavy flavor decays 
was presented. At the time of that work, not all the fixed-order computations were available from which to 
extract the coefficients entering the resummation formula, namely ^3, B2 and D2. A numerical estimate of 
the three-loop coefficient A3 was used, which was obtained in jJU] by fitting the known moments of the 3-loop 
splitting kernels and which has been later confirmed by the exact analytic evaluation |28j . As far as B2 is 
concerned, an approximation based on the q q splitting function at two loops has been assumed, which 
was shown to be rather poor by the subsequent exact computation in . The coefficient D2 was taken 
from its original computation in |32) . There is a misprint in (73(A) in |39) in two terms proportional to Aif32'- 
\og[l — A] — l/21og[l — 2A] has to be multiplied by a factor 2, as found indeed in the recent recomputation of 
the /i-independent terms |38| . With the misprint, the terms proportional to Aif32 would indeed appear at a^, 
while they have to appear only at order a^, as shown correctly in the a expansion of the (73 in eq. (42) of |39| . 

Let us note that the soft terms, i.e. the terms proportional to the coefficients Ai and Di, have the singularity 
closest to the origin in A = 1/2 while the collinear terms, proportional to Bi, have only a singularity in A = 1. 



3.2 Inverse transform to physical space 

The original form factor in u space is recovered by an inverse Mellin transform: 

pC + iOC .AT 

<j(u;a)= ^(l-uy^a^ia), (124) 

where c is a real constant chosen in such a way that all the singularities oi lie to the left of the integration 
contour. The inverse transform can be done to any given logarithmic accuracy in closed analytic form, where 
now the logarithmic accuracy is defined as before but in terms of powers of a and L — log 1 /u instead of 
I ~ logiV. To NNLO accuracy, one can write PHI 

„Lsi(t) + 32(1-) 

r[i-Mr)] ^^' (125) 

where 

r = l3aaL (126) 

and we have defined 

hi{r) = ^[rgiir)] = 9iir) + r g[{T). (127) 
(51] is a NNLO correction factor which can be set equal to one in NLO: 

S^NLO = 1. (128) 

Its NNLO expression reads: 

SJ: = S/S\l^o (129) 

with 



S 



2"5^m|i + (3oag'2{T)tP[l~h,{T)] + ^(3oah[{T){i^' [1 - h^ir)] ~ i^' [1 ~ h,{T)]}\. (130) 



In inhomogeneous terms were included in 6T,, which have been subtracted here. T{x) is the Euler Gamma 
function and 

V'(x) = -^iogr(x) (131) 

ax 

is the digamma function. 



A factor 1 — u ^ 1 has been neglected in our leading twist accuracy. 
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Expanding the r.h.s. of eq. (|125|l up to third order, one obtains the following relations: 



Gi2 

Gn 
G23 

G22 
G21 
G34 

G32 



2 ^' 

-^A2-i/3o(i3i+2Di)-iA?z(2); 
-iB2 + D2) - [Bi + Di) z(2) - Alz{3); 



-^^3 - /3o (B2 + 2 i?2) - y + 2 - Ai A2 z(2) - (5 Si + 7 z(2) 



Ai'iBi +i?i)z(3), 



(132) 
(133) 
(134) 

(135) 
(136) 
(137) 

(138) 
(139) 



where z{a) = 1/"-° is Riemann Zeta function with z(2) = tt^/G = 1.64493 • • • , z(3) = 1.20206 • • • and 

z(4) = 7r''^/90 = 1.08232 • • • . Note that the leading coefficients G23 and G34 involve products of the one-loop 
coefficients Ai and /3o only. The explicit expressions of the Gij read: 



G12 
Gn 
G23 
G22 
G21 
G34 

G33 



G32 — 



_Cf 
27r' 
TGf. 

Cj. f 11 Ca 2n/ 



Cf 
4 7r2 

Gf 
6 7r2 

Gf 
48 7r3 

Gf 

4 7^3 

llz(2) 



G^ I ^ + z(2) ) - 



36 



72 



^ + z(2))+G^ (r^_I_^(2) 
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- 2Gi.z(2) 



905 17 



3z(3) 



48 



Cf ( — + 6z(2) + 3z(3) 



847 Ca^ 77 CAUf 7nf' 



36 



9 



9 



(140) 
(141) 
(142) 
(143) 
;(144) 
(145) 



iZ^ + CAUf — 



108 



Ca Cf + CFnf 



20 z(2) 



27 



Ca' 



1261 llz(2) 



432 



6 



z(2) 



4Gf^z(3) 



(146) 



Cf 

4 7^3 

Ga' 



,/275 _z(2)X 
V 648 n I ^ ^ / 



5399 4z(2) z(3) 
1296 ^ ~ 



6 



+ 



/21893 119z(2) _^ 77z(3) llz(4) 
V 2592 



Cftij I — - 



36 



19 71z(2) 



12 



) +Gf'(-7z(3) + z(4)) + 



48 



36 



+ z(3))+G^Gf (- 



11 685 z(2) 



72 



llz(3) + 5z(4) 



(147) 



The numerical values of the coefficients show a good convergence of the perturbative series also in configuration 
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space: 



Gi2 = - 0.212207; (148) 

Gil = +0.742723; (149) 

G23 = -0.185756 + 0.011258n/ = -0.151982; (150) 

G22 = +0.152206 - O.G12196n/ = +0.115618; (151) 

G21 = +0.628757 - 0.0427065 n/ = +0.500638; (152) 

G34 = -0.189702 + 0.022994 n/ - 0.0006968 = -0.126990; (153) 

G33 = -0.349055 + 0.033368 n/ - 0.0000995 = -0.249846; (154) 

G32 = +0.96117 - 0.09368 + 0.0025974 n? = +0.703506, (155) 



where on the last member of the r.h.s. we have set 71/ = 3. 



4 Distribution in the hadronic variables 



The distribution in the hadronic variables u and w is obtained integrating the triple differential distribution 
l|59|) over the electron energy x = 1 — x. The integration range is 



Xl{w,u) < X < X2{w,u), 



where: 



xi{w, u) 



w u 



and 



X2{w, u) 



W 



(156) 



(157) 



l+u l+u 

Let us use the second method of integration of the triple-differential distribution discussed at the end of sec. 
i.e. let us neglect at first the remainder function. Since the QCD form factor a[u; a(wm},)\ does not depend 
on the electron energy x, the integration only involves the coefficient function: 



I dxC{x, w; a) 



(158) 



We eliminate small terms 0{u) from the integral above by integrating over the range which is the limit w — > 

of (nngi: 

xi{w,0) <x< X2iw,0). (159) 

In fact, these terms 0(u), when multiplied with the plus distributions of u contained in the QCD form factor 
ct(u; a), give at worse terms of the form log u, which miss the 1 /u enhancement and therefore are to be considered 
as "small" . Let us define therefore the coefficient function of the double hadronic distribution as: 



Ch {w; a) — dxC {w, x; a) , 
Jo 

having the usual a expansion 

Ch{w; a) = cPiw) + aG|^^^(w) + dj}\w) + 0{a^). 

One easily obtains: 

C'Piw) = 2w2(3-2w); 

C'j^\w) = —w'^\-i9-4w)logw + 2i3~2w) 



(160) 



(161) 



(162) 



35 

Li2 {w) + log w log(l — w) — 



-"^^The latter are actually the integration regions for e+e — » <7 + <? + 9 with massless quarks. 
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The first two orders of tlie coefficient function vanish as uP' for w — > 0, implying a suppression of the states 
with a small hadronic energy (i.e. with a small hard scale), as anticipated in the introduction. 

The resummed distribution in the hadronic variables u and w then reads: 
1 fPV 



= Ch[w\ a(™fc)] o- [w; a{'wrrn,)\ + dulu^w] Q!(mb)] , (163) 



r dudw 

where the remainder function has an expansion analogous to the one in the triple differential distribution: 

dH{u,w;a) = Q;d^^(u,w) + d^^\u,w) + O(a^). (164) 

Expanding to first order the above distribution and comparing (matching) with the known 0{a) distribution, 
the following remainder function is obtained — an over-all factor Cp /t^ is omitted: 

Ji), , 4?«^Iogu 4w^ (3 - 2 w) log u 32 - 10 log u 

dW(w,U)= - =- ^ ; 7j h 

3 (14 - 6 w - 5 - 2 (3 - 4 w) log M 
^ (l + «f ^ 



20 w3 (2 + w) (1 - 2 w) - w-* (9 - 18 w - 2 w^) log u 
6Aw^ + 2w^ (3-6w-4w2) log?/ 



+ 



{l + uf 

4 (10 - 15 w - 2 w2) -w^ (12 - 13 If; - 6 log 
(l + «f 



(165) 



Eq. (|163|l provides a complete NLO rcsummation of the distribution in the two hadronic variables u and w, 
from which the distribution in any other pair of hadronic variables can be obtained by a change of variables. 
One can insert in eq. H163|) the NNLO form factor a, whose properties have been discussed in sec. ©, allowing 
an approximate NNLO resummation. In fact, for a complete NNLO resummation, one also needs the second 

(2) (2) 

order corrections to the coefficient function C}j {w) and the remainder function d\j {u,w), which are unknown 
at present. 



5 Hadron energy spectrum 

The distribution in the total hadron energy w is obtained by integrating the distribution in the hadronic variables 
(|163|) . The integration range in u is: 

max(0,w-l) < u < 1. (166) 

Since the coefficient function Ch{w; a) does not depend on u, the integration only involves the QCD form 
factor and the remainder function: 

— — — = Ch (w; a) |l — 6{w — 1)1] [w — 1; a(wmf,)] | + / dw dH{u,w; a), (167) 

T dw I J Jmax(0,tu-1) 

where S(u; a) is the partially-integrated form factor defined in section 



Because of the 9{w — 1) multiplying T.(w — 1; a), there are large logarithms only for it; > 1, as anticipated in 
the qualitative discussion in the introduction. We may therefore consider the parts of the spectrum for w < 1 
and w > 1 as two different spectra, merging in the point w — 1. Let us consider the simpler case w < 1 first. 
Since, as already noted, there are no large logarithms, no resummation is required and the 0{a) fixed-order 
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result coincides with the NLO one. There is no QCD form factor and therefore there is no way to distinguish 
between the coefficient function and the remainder function. The spectrum for w; < 1 can then be written as 
an ordinary a expansion: 



1 dT 

2f = ^("'^ 



[w < 1), 



where: 



L{w; a) = L^"\w) + a L^^\w) + a' L^^\w) + 0(a3). 
The first two orders read [7| |S] : 



(168) 
(169) 

(170) 
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— +Li2(l - w) 



720 



(4w'' -42u;3 + 585^2 -3720w; + 4860 + 1440?/; log w- 3240 log w) [. (171) 



Let us now consider the more interesting case w > 1, where resummation is effective and one has to keep 
the resummed form of the distribution in l|167l) . In a minimal scheme we have to subtract small terms from 
the first term on the r.h.s. of eq. H167|l . since the form factor must contain large logarithms only. This is done 
setting w = 1 in the argument of the coupling entering the form factor S as well as in the coefficient function 
Cjy, obtaining the simpler expression: 

= (1; «) {l - S [ii; - 1; a(mf,)] } + •••, (172) 

where the dots denote terms not containing large logs of w — 1. Let us prove the legitimacy of the transformation 
from (|167|) to (|172|l . As far as the argument of the coupling is concerned, we expand the QCD form factor S in 
powers of aiwurib). One obtains terms of the form 

a(wr7ifc) log^(w— 1) = Q:(mf,) log^(iy — 1) — 2/3oa(mb)^ logwlog^(u' — 1) + 4/3Qa(mb)'^ log^ wlog^(w — 1) + • • • , 

(173) 

where on the r.h.s. an expansion of a{wm}j) around the point w = 1 has been performed. All the terms on 
the r.h.s except the first one vanish for w ^ 1^, therefore they are not large logarithms and can be dropped. 
The only large logarithm is the first term on the r.h.s., which is obtained by setting w = 1 vti the coupling in 
the original expression on the l.h.s. All this implies that the coupling can be evaluated in the infrared-singular 
point 1/7 = 1. As far as the coefficient function is concerned, one just notices that the neglected terms, 

Cniw- a) - Ch(1; a) ] {l - S [u. - 1; a(mb)] }, (174) 

are again vanishing for w — > 1"*", because Cr (w; a) — Ch (1; ct) = 0{w — 1) and therefore can be neglected in 
this limit. 

In NLO one has also to add a remainder function to be determined via a matching procedure. That, as 
already discussed in other cases, is in order to take into account also the region w — 1 0{1). One then has 
the resummed expression: 

^ Cwia)[l - ^[w-l;aimb)]} + Hiw;a) [w > I), (175) 

where we have defined: 

Cwia) = Ich{1; as). (176) 
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The coefficient function and the remainder function have a standard a expansion: 



Cwia) 
H{w; a) 



= 1 



- 0{a^). 



(177) 
(178) 



The first order correction to the coefficient function reads 



a 



(1) 

w 



1.15868. 



(179) 



Note that is negative and has a rather large size; for a(mf,) — 0.22 it gives a negative correction of 

« — 25%. By using the matching procedure described at the end of sec. (O, we obtain: 

i7(i)(w) = — ( - i (2w + 1) (u; - l)2log2(i(; - 1) - i(w - l)(2w2 - w -4)log(w - 1) + 
TT L 2 3 



2Li2 (1 - w) + 21og(w - l)logw H 

6 



+ ^(2 - w){4w'^ - Mw^ + 517w^ - 2M6w^ + 3798w + 1248) |. 



(180) 



The above function is positive in aU the kinematical range 1 < w < 2 and goes to zero for w 
on the basis of the vanishing of the phase space in this point. 



2, as expected 



Let us make a few remarks about eq. (|175f) . If we expand the r.h.s. of cq. (|175|l in powers of a, we find that 
c|^'' only appears in order — this occurs because the form factor multiplying the coefficient function is in 
this case 1 — S = 0{a) and not S = 0(1). At present, only a full 0{a) computation is available, implying that 
Cw'' cannot be determined by the matching: only the remainder function can be fixed by this procedure. The 



(1) 



value of 

is however another method to fix C^: we require that the resummed spectrum is continuous in w 



E(u;- 



^ came out "automatically" as a consequence of our resummation formula (see. eq. (|176|l V There 

1. Since 

1) ^ for w 1+, we obtain the equation: 



1 + aLW(l) = 1 + a\clP + H^'\^) 



(181) 



to be solved in 



a 



(1) _ ^(1) 



- ^ >il) - H^^'il) (182) 

and giving again the value H179|l . The condition of continuity of the resummed spectrum in w = 1 is very 
reasonable from the physical viewpoint and it is remarkable that the two methods give the same value for the 
coefficient function. 

Even though we are considering a differential spectrum, its resummation involves, as we have explicitly seen, 
the partially integrated form factor. S usually enters event fractions in expressions of the form 



i?(y; a) = C(a)S(?/; a) + D{y; a), 



(183) 



with a remainder function vanishing for y —^ 0, where y is a general kinematical variable entering the large 
logarithms log 1/y. In the case of the hadron energy spectrum, its resummation is different from (|183l) because 
it involves the combination 1 — E instead of E: there is an additive constant, namely one, which makes the 



spectrum non vanishing for w V 
remainder function H(w]a) for w - 
this condition, because: 



as it should. It seem however reasonable to impose the vanishing of the 
1+ also in this case. The previous factorization scheme does not satisfy 



iJ«(l) 



Cf_ / 2587 



TT 

~6 



(184) 
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We can construct an improved scheme satisfying this condition by introducing two coefficient functions instead 
of one: 

— — Cwi M {l — Cw2 (a) S [w — 1; a(TOfc)] + Hiw; a)\ (improved scheme, w > f), (185) 

2r dw L J 

where the new remainder function, vanishing in w = f , reads: 

H{w] a) = H{w; a) - ff(l; a). (186) 
The coefficient functions have the usual fixed-order expansions: 

Cwi {a) = 1 + a&^\ + a^&^\ + ©(a^); (187) 
Cw2{a) = 1 + aCl^l, + a"c|^^2 + 0(a3). (188) 

By imposing the continuity in w = 1 as in the previous scheme, we obtain for the first coefficient function at 
first order in a: 

CW=L«(1) = -^1^^-0.331868. (189) 
The second coefficient function is obtained by imposing the usual matching with the first order computation: 

c'^\ - H^'\i) = - V ( W - y) - - ^''^^ 

The improved resummed expression (|185ll is positive in all the kinematical range 1 < w < 2 and vanishes for 

We can compare the hadron energy spectrum for w > 1 given in eq. (|175|l or in eq. (|185|1 with the hadron 
mass distribution in the radiative decay Q given in eq. (|50|l . The hadron energy distribution contains S, i.e. 
just the integral of the form factor a entering the radiative decay spectrum. The hadron energy spectrum 
is therefore a very good quantity on the theoretical side — it is exceptional in this respect — being directly 
connected, via integration, to the radiative decay. By that we mean that the connection between the two 
spectra only involves short-distance coefficients. As show in 5 , this is to be contrasted with the case of other 
single-differential spectra. 

5.1 Average energy 

As discussed in the introduction, the infrared singularity in w = 1 of the 0{a) spectrum is integrable, so one 
can calculate directly the average hadronic energy as a truncated expansion in a: 



Iw) = — 
^ ' 10 



^ aCp 137 
^ TT 840 



0.71. (191) 



The 0{a) correction is very small, of the order of 1%, due to a large cancellation between the contribution for 
w <1, which is negative, and the one for w > \, which is positive. Setting for instance m;, ^ ms one obtains 
in leading order: 

(i^x) = = 1.843 GeV (192) 

with a tiny first-order correction of -I- 26 MeV. This quantity can be directly compared with the experimental 
value. In the radiative decay ((SJ there is a larger final hadronic energy: in lowest order 

{Ex)B^x,y = ^ = 2.634 GeV. (193) 

The average hadronic energy is ^ 30% larger in the radiative decay than in the semileptonic decay, in line with 
the qualitative discussion about the differences of the two decays given in the introduction. 
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5.2 Upper cut on hadron masses 



In experimental analysis an upper cut on invariant masses 

mx < mx (194) 
is imposed in order to kill the large background from semileptonic b c transitions. Let us define: 



(195) 



In practice, mx — 1.6 -f^ 1.8 GeV, so we can assume fc < 1. A leading order evaluation of the spectrum with the 
above cut gives: 



1 dT 
where 



w2(3 - 2w) { d{k -w)+ d{w - fc)S 




l-y/l-(k/wy\ 
l+^l-(fc/to)2 ■ 



a(wmb) 



WM 



9{w — 1)S [w — 1; a{mb)] ^ w < Wm 

w > Wm 
(196) 



(197) 



is the maximal hadronic energy above which the spectrum vanishes; as expected on physical ground, cutting 
large hadron masses also acts as an upper cut on hadron energies. The spectrum is continuous in w = wm and 
it develops large logarithms for fc — > 0. Let us observe that the argument of the first QCD form factor E has a 
similar form to the variable u defined in eq. (|36|l . In fact. 



mx 
Ex 



(198) 



is the analogue of the variable 4y with y defined in eq. (|57|l . 



6 Distribution in hadron and electron energies 

In this section we derive the distribution in the hadron and electron energies w and x by integrating the triple 
differential distribution l|59|) over u. In general, there are two independent energies in the semileptonic decay 
lO. That is because the hadronic final state is basically a pseudoparticle, i.e. a single entity possessing an 
energy Ex and a (variable) mass mx- We have therefore 3 particles/pseudoparticles in the final state and 3 
energies, related by energy conservation: 

Xe + x^ + w ^ 2, (199) 

where 

x^ = — - (200) 

mb 

and we have written Xe instead of x for aesthetical reasons. Since the neutrino energy is not usually measured, 
let us take as independent energies the electron and the hadron energies. We have to integrate over u in the 
range 

' w 



max[0, w — 1] < u < 



(201) 



As in the previous section, let us use the second method of integration, i.e. let us omit at first the remainder 
function. Since the coefficient function C{x,w] a) does not depend on u, the integration only involves the form 
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factor a — a complementary situation with respect to the one in the previous section — and we obtain: 



1 d^Y 
r dxdw 



C'[x, w; a(mf,)]| 6{2x — w)T, 
- 9(w-l)T. 



W — X 

— — — ; a{w mt) 

X 



+ d{w - 2x)E 



w — X 



— ; a{w mjj) 



+ 



w ~ 1: a 



} 



(202) 



where the dots denote non logarithmic terms to be included later. The decay ((TJ involves an hadronic subprocess 
with a heavy quark decaying into a light quark evolving later into a jet. Hadron dynamics is therefore symmetric 
under the exchange of the electron and the neutrino momenta, since it is "blind" to W decay. That is clearly 
seen by expressing w through x^, by means of eq. (|199|l : 



1 d^T 



r dXedXy 




; a((2 - Xe - Xy)mb) 
-; a((2 - Xe - Xy)mb) 



l-x 

(1 — Xe ~ x^)Yj 1 ~ Xe — Xy] a{mi,) '> + 



(203) 



Soft-gluon dynamics — i.e. the expression above in curly brackets — is symmetric under exchange of Xe with x,^. 
The coefficient function C[xe,x,y; a{mb)] however is not symmetric under the exchange of the lepton energies 
because it does depend on the whole process, involving the decay of the W boson into the lepton pair, and not 
only on the hadronic subprocess. 

To proceed with resummation, however, let us go back to the more familiar variable w, i.e. to eq. I|2U2I) . 
Large logarithms can in principle be obtained by sending to zero the argument of any of the QCD form factors 
E's entering 1)2021) . i.e. in the following three cases: 



1. w 



0; 



2. X 



0; 



3. w 



1" 



(204) 



The coefficient function C\x, w\ a(m(,)] vanishes in the first limit as 0{w — a;), implying that in this case there 
are actually no large logarithms. This limit corresponds to —f mb/2, a point where the tree-level spectrum 
vanishes suppressing soft-gluon effects. The only relevant limits are therefore the second and the third ones. It 
is therefore natural to write a factorization formula dropping the form factor not associated to large logarithms: 



1 d^r 



= c 



[x,w] a(mb)] |l] [x/w; a{wmb)] — 9{w — 1) S w — 1; a(rnb) | + 



r dxdw 

We have taken the limit a; — > in the theta functions containing x in the argument. 



(205) 



Let us consider separately the cases w < 1 and w > 1. In the simpler case w < 1 there is a single form 
factor and one can write a factorized expression of the form: 



1 d^T 



— Cl {x, w; a) S [x/w; a{w rrn,)] + d^{w, x; a) 



T dxdw 

We require that the remainder function vanishes for a; — + 0: 

]mi d<^{w,x; a) = 0. 



(w < 1). 



(206) 



(207) 



^^Notc that this case is a "comphcation" of the analogous case for the single distribution in w, where the integration over x has 
been made and therefore there are no large logarithms of x. 
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The coefficient function {x,w; a) can be taken as: 

Cf\w,x) = I2{w~x){l+x-w); (208) 



C^l\w,x) = —I2{w-x){l + x~w) 



T ■ / \ 1 1/1 \ 3 , wlogw 35 

Li2(u;) +logwlog(l - w) - -logw - _ - — 



(209) 



In cf \'w, x) we have put the factor 12{w — x){l + x — w), equal to the spectrum in lowest order, in order to 
have a vanishing remainder function in 0(q;°): this is a non minimal choice, since the minimal choice would 
imply to set a; = in the coefficient fimction. We have inserted a similar factor also in C^j^\w,x), in order to 
have a simple multiplicative form of the correction As in previous cases, by matching with the full 0(a) 
result m, we determine the remainder function 

d<{w,x; a) = ad^^\w,x)+ d^^\w,x) + 0{a^). (210) 
Omitting the over-all factor Cf/tt, we obtain for the leading contribution: 

d^^\w,x) = -^(w-^) x(- 210 + 280w-10w2 + 2w^-60x-125w^-7w^x + 15^^ + 
+ 32wx2 „ 15^) 

+ ——^ ( - 45 1« + 60 u;^ - 20 + 10 w'' - 6 + - 15 X + 135 w X - 255 ty^ X + 

5 (-l + w)^ 

+ 85w^x + 25u;^x - bw^x- 15x^ + 45 wx^ + 75w^x^ - 85w^^^ + IQw'^x'^) \ogw + 

- 6 (—1 + w — x) (w — x) log' w + Q (—1 + w — x) [w — x) log^(w — x) + 

- -{w -x) (45 - 15 w + 5 - 5 + + 15x - lOwx + 15w^x - 4w^x + 5^^ + 
5 ^ 

- Ibwx'^ + 6u;^ + 5x^ - 4u;x^ +x^) \og{w -x) + 

- -x{m-l%Qw + l2Qw'^ + mx-lbwx-^bw^x + ?>x'^ -2Qwx'^ + lQw'^x'^ + 
5 ^ 

+ bx^ -bwx^ + x^) logx + 

+ 12 [—1 + w — x) {w — x) logw logx — 12 {—\ + w — x) {w — x) \og{w — x) \ogx . (211) 

Let us now consider the case w > 1: 
1 d^T 



r dxdw 
where we have defined 



= C\ 



{x, w; a) |l] [x/w; a{w mb)] — S Aw; a(mf,) | + • • • (w > 1), (212) 



Aw = w - 1 > 0. (213) 



There are two form factors and large logarithms can be obtained in the following three kinematical configura- 
tions: 

1. a; ^ Aw ~ 1: large logarithms of the form a" log*^ x have to be resummed; 

2. Aw ^x-^1: large logarithms of the form log Aw ~ logx have to be resummed; 

3. Aw <C X ~ 1: large logarithms of the form log Aw have to be resummed. 



■^"We could have taken as coefficient function the original one C{x,w, a) as well, which however does not always contain the 
factor 12{w — x)(l -|- x — to). 



29 



The first case is kinematically forbidden because 

Aw; < X. (214) 

The second case does not give large logarithms because the coefBcient function C{x, w; a) vanishes linearly in 
this limit: 

C{Xx,l + XAw; a) ^ 0{X) for A ^ 0. (215) 

The only relevant limit is therefore the third one, implying that one can drop the form factor T,{x/w; a). We 
propose then a resummed form for this distribution which is a generalization of that one for the hadron energy 
spectrum: 



1 £T 



T dxdw ^^^1 ~ CXW2 (x; a)S Aw; a{mb) | + d>{Aw,x; a) {w > 1). (216) 

We require that the remainder function vanishes for Aw; 0+: 

lim d>{Aw,x; a) = 0. (217) 

Aiu— >0+ 

The first coefficient function is obtained by imposing the continuity of the spectrum for w; — > 1 from both sides 
w; < 1 and w; > 1 and for any x . We obtain: 

Cxwi{x; a) = Cl{x,1; a)T,{x; a) + d<(l,S; a). (218) 

We can expand in the above equation S {x; a) in powers of a (up to first order) because the coefficient function 
for w = 1, Cl (^, 1; ol), vanishes linearly for x — > 0, killing the large logarithms in the form factor. We then 
obtain: 

Cxwi{x; a) = c'^Uix) + aC%,ix) + cf^^{x) + ©(a^) (219) 



where 



= l2{\-x)x; (220) 
Cwi(^) = ^ (-587+192^-47^2^ 15 x^) - \x (105- 105^-5^^+^) \ogx + 

- ^(1-^) (31 + 16x-4^2 +x'') log(l-x)-6 (1-x) X log2(l-x) + 





+ 12 (1 -S) S log(l -x) logs- 6 (1 - x) X log2s+ 12 (1 - x) ^0(2) [>. (221) 



The second coefficient function Cxwi{x; a) is obtained by matching with the fixed-order distribution in the 
limit Aw; ^ 0+: 

Cxw2{x; a) = C%^{x) + aC%^{x) + C%^{x) + ©(a^) (222) 



with 



Cw2(^) = 1; (223) 

Cxwo{x) = — I— (62- 192^ + 47x^-15^^) + — (31 + 16x-4x2+x3+x*) log(l-S) + 
XW2\ ^ TT 1^ 120 ^ ' 60S ^ / &v / 

+ - \og^{l-x) + , ^ _, (105 - 105S- 5^2 +x'^) logs - log(l-S) logS + 

2 60 (1 — a;) 



+ \Wx\. (224) 



•^'^This continuity condition, which involves a single point w = 1 for the hadron energy spectrum, involves in this more complicated 

case the line {w = 1, x). 
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The remainder function dy{w,x] a) is obtained by matching with the fixed-order distribution for Aw ^ 0(1): 

dy{w,x; a) = ad^yHw,x) + d^^\w,x) + O(a^), (225) 
where, omitting the overall factor Cf/tt: 

d^^\w,x) = ^(-1 + u;) (1-x) (- 75 + 142u)-7w^ + 2w^-212x-105iyx-9u;^x + 

+ 132x2+39wx2-47x^) + \{-l + w) ( - 4 - 99 + - 4 + + 140x + 

5 

+ \2Qwx + lb'w^x~bur"x-Qbx'^~ 65 wx^ + lOw^x^) log(-l + w) + 

- 6 i-l + w) (w-2x) log^(-l + u;) + i(l-S) (Sl + Wx-Ax^+x^ + x^) log(l-x) + 

5 

+ 6 {1 — x) X log^(l — x) + 6 (—1 + w — x) {w — x) log'^(w — x) + 

- (-1 + w) {-'2lw + 30x + 2Awx-12x^ -9wx'^ -2x^ + 2wx^ -x^) logx + 

- 12 (1 — x) a; log(l — x) logx + 6 (—1 + w) {w — 2x) log^ x + 

- ^{w -x) (45 - 15 w + 5iu^ - 5iy^ + + 15x - 10u;x + 15w^x - 4w^x + 5x^ + 

- I5wx'^ + 6 + 5x^ - Awx^ +x^) \og{w -x) + 

- 12{—l + w — x){w — x)logx\og{w — x). (226) 

To summarize, we have presented a complete NLO resummation of the distribution in the hadron and electron 
energies w and x, which is a generalization of the resummation of the hadron energy spectrum of the previous 
section. Resummation takes a different form in the cases w < 1 and w > 1. In the first case there is a series of 
threshold logarithms of the form 

a" log'^ - {w< 1), (227) 

w 

while in the second case the infrared logarithms are of the form 

a"log'=(w-l) (w>l). (228) 

Unlike the distribution in sec. 01 we have here a differential distribution involving the partially-integrated form 
factor S. 



7 Conclusions 



It is a rather old idea that semi-inclusive B decays can be related to each other because of some universal 
long-distance component |2S]. We have presented in this paper a critical analysis of this idea, based on a 
resummation formula for the triple differential distribution in the semileptonic decay Q . Long-distance effects 
manifest themselves in perturbation theory in the form of series of large infrared logarithms, coming from the 
multiple emission of soft and/or coUinear gluons. The universality of long-distance effects has therefore to 
show up in perturbation theory in the form of identical series of large logarithms in different distributions. 
Semi-inclusive B decays have been defined in all generality as decays of the form 

B ^ Xq + (non QCD partons), (229) 

in the kinematical region close to the threshold mx = 0, i.e. for 

mx <t: Ex. (230) 
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We have shown that semileptonic distributions are naturally divided into two classes. 

The first class contains distributions which are not integrated over the hadronic energy Ex and consequently 
have a long-distance structure similar to the one in radiative decays (jSJ. These are the (simpler) distributions 
to attack and have been treated in this paper. We have resummed to next-to-leading order: 



1. the distribution in the hadronic energy Ex and in the variable u defined in sec. (O, which is basically the 
ratio in^xHAE"^), i.e. the hadron invariant mass squared in unit of the hard scale; 

2. the hadron energy distribution, which is a case of the so-called Sudakov shoulder. This is the only 
single distribution which can be related to the radiative decay via short-distance factors only. The large 
logarithms which appear in this distribution are indeed equal to the ones which appear in the radiative 
decay (jSj). We have studied in detail the relation between the hadron energy spectrum and the photon 
spectrum in the radiative decay. It is remarkable that the large logarithms in the hadron energy spectrum 
occur at Ex — '77-b/2, i.e. when the hard scale Q — 2Ex equals irib, as in the radiative decays; 

3. the distribution in the hadron and in the charged lepton energies, which contains two different classes of 
large logarithms according to the cases w < 1 or w > 1. The resummation of this distribution is the most 
complicated and is a generalization of the resummation of the hadron energy spectrum. 



The second class contains semileptonic distributions in which the hadronic energy is integrated over, such as 
for example the hadronic mass distribution or the charged lepton energy distribution. These distributions have 
a complicated logarithmic structure, which is not simply related to the one in the radiative decay and there is 
not a pure short-distance relation with the radiative decay spectrum. The resummation of these distributions 
to NLO is presented in 0. 
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